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1 Introduction 

CN 

i— | It is well known that for an arbitrary convex solid in three dimensions, the aver- 

age projected area is one-fourth the surface area. Perhaps less well-known is the 
rather touching story of how this theorem originated. Karl Schwarzschild, a Ger- 
man astronomer primarily famous for his black hole solution to Einstein's gen- 
eral relativity, proved the theorem near the turn of the twentieth century. How- 
ever, perhaps due to his untimely death at forty- two in 1916, the theorem was 
£h not well-known. So it remained for his son, astronomer Martin Schwarzschild, 

l — 1 noted for his contributions to understanding stellar evolution, especially red gi- 

ant stars, prove the theorem after the second World War, unaware of is father's 
work (Knapp, personal communication). One can only imagine his surprise to 
learn from a colleague that his father had already proved it decades earlier! In 

^\ the present work, I generalize this theorem to higher dimensions. 

What is meant by average projected area? There are two equivalent, conve- 
nient ways to imagine this. The first is, hold the solid's position fixed, and shine 

ON a light beam behind it toward an observer in front of it. Vary the light beam's 

direction (and consequently the observer's location, since two points define a 
line and so the observer must move to see the shadow) through all possibilities, 
• • measure the area of the shadow cast at each, and average. The second is, hold 

the light beam (and observer's) position fixed and rotate the solid through all 
possible orientations, keeping its spatial location fixed. Measure the shadows 
?_i and average. Either method yields the same result - as one would expect, since 

they correspond to the same process but observed either in the rest frame of 
the solid or of the observer. 

So we have set out what average projected area is. But why should one care 
about the average projected area of a convex solid? Applications are many. In 
astrophysics, dust is along the path light takes from distant stars or galaxies 
to us, and so extinguishes some of the light. The dust grains are irregularly 
shaped and oriented, so it is necessary to correct for their average extinction of 
the incoming light. This theorem provides the most convenient way to do so. 
The next step is complex numerical simulations; see e.g. Draine 2011. 
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The theorem also readily allows one to calculate the temperature of dust 
grains (or asteroids) illumined by a star (assuming the grains or asteroids are 
convex). The energy received by the dust grain is just proportional to its pro- 
jected area (projected onto the light source): clearly, it cannot receive energy 
from any part of its surface not "visible" to the illumining star. However, the 
particle can radiate energy away proportionally to its surface area. Setting en- 
ergy absorbed per unit time equal to energy radiated per unit time (i.e. the 
luminosities L a b s and L ra d) determines the equilibrium temperature: 

irad = AsfSB^grain = £ a bs = (A pI0 j) ^SB^* ^ 1 (1 - a) . (1) 

As is the surface area of the grain, ctsb the Stefan-Boltzmann constant, T gra ; n 
the grain temperature, (A pl0 j) the average projected area of the grain, the 
star's temperature, i?* the star's radius, d the distance between the grain and 
the star, and a the albedo, a measure of how reflective of light the grain is 
(the higher a, the less energy the grain absorbs). Clearly, without the ratio 
(^4proj) /As, it would be impossible to calculate the temperature of the grain. 

The average projected area theorem also has more down-to-earth applica- 
tions. It is a starting point for determining the 3-d shape of micro-particles 
using a 2-d map of their projected areas (see Vickers & Brown 2001, Brown et 
al. 2005). The theorem has military applications as well: Saucier (2000) uses 
it in discussion of simulating debris fragments behind armor and characterizing 
the penetration potential of projectiles. Further, it is a useful approximation 
when considering particle transport in 3-d; see Glassner 1995. The theorem also 
appears as a step in estimating the sphericity of fruits and of polymer beads in 
gels; see e.g. Houston 1957 and Nussinovitch 2010. Finally, the theorem has 
applications in ray tracing for computer graphics (Hanrahan 2000). 

The main purpose of this paper is to generalize this theorem to higher- 
dimensional convex solids. I show how to compute (and prove) the analogous 
factor to one-fourth (in (^4 pro j) = \As : with As the surface area and angle 
brackets denoting average) for an arbitrary dimension: for example, one can 
use my results to compute the 7-d average projected area of an arbitrary 8-d 
convex solid. I also touch on the limit of this factor as the dimension becomes 
infinite. 

The paper is structured as follows. For completeness' sake, I begin by prov- 
ing the theorem in 3-d; my treatment is original to this paper, but proofs can be 
found elsewhere, e.g. in Hildebrand 1942. In Section 3, I prove a method of cal- 
culating the analogous ratio for higher dimensions that uses hyperhemispheres. 
Section 4 applies this method to compute the ratio of average projected area 
to surface area as a function of the dimension of the space, provides a table of 
values, and examines the limit of the ratio as the dimension becomes infinite. I 
continue in Section 5 by presenting a simple, memorable recursion relation for 
this ratio between dimensions d and d+1. In Section 6, I derive from the recur- 
sion relation an explicit formula for the ratio in terms of elementary operations 
(the formula I supply in Section 4 is in terms of gamma functions). 
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2 Proof in 3-d 



Here we prove that for a convex solid in three dimensions, the average projected 
area is one-fourth the surface area. 

Consider an arbitrary convex solid, and imagine an infinitesimal chunk of 
surface area dA. Let dA be the unit vector normal to the surface. Now imagine 
surrounding the solid with a sphere, and that we, the observer, sit at a point 
on the surface of that sphere. If the vector dA points such that it intersects 
the hemisphere whose apex is our location, then we will see it. Since this 
hemisphere's position is uncorrelated with the direction of dA for an arbitrary 
chunk of surface area dA, on average dA will point in the right direction for us 
to see it half the time. 

We now restrict our attention to the case where dA does point into "our" 
hemisphere. Consider its projection along the line of sight. Denote the angle 
between dA and the line of sight by 9, and realize that because we have restricted 
ourselves to a hemisphere, 9 goes from zero to tt/2. 

We wish to integrate cos 9 to find its average value, as cos 9 gives the projec- 
tion of dA onto the line of sight. When so doing, we will use as our Jacobian for 
spherical coordinates sin9d9. Using this Jacobian builds in a certain convention 
about how 9 is defined: namely, that it goes from zero to tt. (Basically, think 
of the line of sight as like the z-axis in spherical coordinates, so 9 is the usual 
angle between the z-axis and dA.) 

And we must ignore the contribution of 9 between tt/2 and tt because that 
would correspond to dA pointing into the hemisphere that is not visible to us. 
So 9 runs from zero to tt/2. 

We now have 

/ - \ r' 2 1 

(dA-zj= J cos (9 sin 6W6> = -. (2) 

(Notice that this integral is really just the differential projection of the solid's 
face onto the hemisphere we are sitting at the apex of, normalized (divided by) 
the surface area of the hemisphere.) 

Here, the picture is that the orientation of the solid's face is changing while 
the observer's position is held fixed. In this picture, we are working in the 
rest frame of the observer. As we noted in the Introduction, it is equivalent to 
consider the solid's face fixed and vary the observer's position with respect to 
it - in other words, to work in the rest frame of the solid. The only quantity 
of physical significance is the relative angle between the line of sight and the 
oriented unit area vector dA. 

Taking into account the additional factor of 1/2 because dA only points into 
the hemisphere visible to us half the time, we have 

(A proj , 3 - d ) =J(dA-z)dS= \a s , (3) 

where S denotes an integral over the surface of the solid and As is the solid's 
surface area. So this is the result for three spatial dimensions. 
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Before we move on, note that, although we derived equation ([2]) in the rest 
frame of the observer, working in the rest frame of the solid will prove more 
fruitful for generalizing the theorem to higher dimensions. For notice that, in 
this picture, as the observer's position varies (at fixed distance from the solid's 
face), it describes a hemisphere. So we simply wish to find the normalized 
(divided by the surface area of the hemisphere) projection of the solid's face onto 
this hemisphere. But the projection operator, or dot product, is commutative, 
so this is just the projection of the hemisphere onto the solid's face - which is 
simply a disc bounded by the circle created where the hemisphere intersects the 
plane of the solid's face. 

As noted, to normalize, we divide by the surface area of the hemisphere. We 
also multiply by the factor of 1/2 because only half the time is the face oriented 
into our hemisphere - or, equivalently, we could normalize by dividing by the 
surface area of the sphere, since that is twice that of the hemisphere. 

3 Proof of method for higher dimensions 

As we pointed out in the previous section, for the 3-d case, considering the 
rest frame of the solid's face and projecting onto the hemisphere defined by the 
different possible orientations of an observer at fixed distance is a fruitful way to 
calculate the average projected area. As we also pointed out, this projection can 
be easily evaluated by finding the normalized projection of the hemisphere onto 
the solid's face. Both of these ideas are also valid in arbitrary higher dimensions. 

In higher dimensions, we simply consider the hyperhemisphere generated by 
varying the observer's orientation (at fixed distance) with respect to the solid, 
and look for the normalized projection of the solid's hypersurface onto it. This 
projection is just equal to the normalized projection of the hyperhemisphere 
onto the hypersurface, which bounds it (imagine the hemisphere as an upside- 
down bowl and the solid's face as the table it's on). And the projection of a 
hyperhemisphere of dimension d onto the bounding hypersurface in d dimensions 
is simply the hypersphere of dimension d— 1. So we see that 

(<M proj)d ) = ^= L , (4) 

where Vd-i is the volume of the unit hypersphere of dimension d — 1 and 
S*h, d its surface area (the factor of 1/2 is because the solid is only pointing into 
our hemisphere half the time). To continue the upside-down bowl analogy, here 
Vd~i is the area of the disc defined by the bowl's edge on the tabletop, and 5h. d 
is just the normalization factor. 

4 Calculation in higher dimensions 

Here, we use the observation of the previous section to derive the analog of 
equation (13J) for dimensions other than 3. As equation Q shows, the differential 
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average projected area of a d- dimensional convex solid is just the volume of the 
unit d— 1-sphere normalized by double the surface area of the unit d-hemisphere. 
Hence we seek the factor k(d) relating (A plo - h d) to As, the surface area of the 
convex solid. As we will see below, it is just Vd-i/S^, where V^-i is the volume 
of the unit d — 1-sphere and Sd is the surface area of the unit d-sphere. We have 

(A proj , d ) = / (dA proj , d ) dS = = k(d)A s , (5) 

JS b d 

where the second equality is from using equation Q for (cM pro j, d) and 
replacing 25h, d with Sd- We can check that the result is correct for 3-d: with 
d = 3, V d -i = V 2 =tt and S d = 4tt. 

Now. for a sphere, 

Vi-i = (6) 



and 



27r d/2-l/2 1 



r(fd-§) MdV* 



(Hypersphere, Wolfram), where we have defined 



We can thus see that 



k(d)= r( ,\, Af ■ (9) 
y/Ti{d - l)M d 

Equation ^ is easily evaluated by e.g. Mathematica; we provide a table of 
the first thirty-two values, as well as a plot illustrating the fall-off in k{d) with 
increasing d. 

It is possible that measurements of this factor could tell us something about 
the number of dimensions we live in. String theory suggests we may live in 
9 (supersymmetric theories), 10 (M-theory), or 25 spatial dimensions (bosonic 
string theory); for reviews, see Kiritsis 1998, Polchinski 1998, or Schwarz 2000. 
However, in string theory these extra dimensions are often constrained to be 
on very small scales - for instance, in Kaluza-Klein theory, the extra dimension 
must be smaller than 10 -16 cm (Yagi et al. 2011), and in Randall-Sundrum 
models (4 spatial dimensions, see Randall & Sundrum 1999 I and II), constraints 
are on the order of micro (10 -6 ) meters (Kapner et al. 2007). It is therefore 
not clear, and conceivably a topic of further research, whether extra dimensions 
would have effects measurable through a value of k(d) different from fc(3) = 1/4. 

Finally, it is interesting to consider the limiting behavior of k(d) as d tends 
to infinity. As might be intuited from the plot, k{d) — > as d — > oo. To prove 
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Average projected area/surface area 
as a function of dimension 
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Figure 1: This plot shows that k(d) falls off as d increases; to leading order 
k(d) oc 1/y/d, as equation (10 1 indicates. 
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Table 1: Ratio of average projected area to surface area (k(d)) as a function of 
dimension d 



d 


k{d) 


d 


k(d) 


2 


.318 


18 


.095 


3 


.250 


19 


.093 


4 


.212 


20 


.090 


5 


.1875 


21 


.088 


6 


.170 


22 


.086 


7 


.156 


23 


.084 


8 


.146 


24 


.082 


9 


.137 


25 


.081 


10 


.129 


26 


.079 


11 


.123 


27 


.077 


12 


.118 


28 


.076 


13 


.113 


29 


.075 


14 


.109 


30 


.073 


15 


.105 


31 


.072 


16 


.101 


32 


.071 


17 


.098 


33 


.070 



this, we consider a series for k(d) about d = oo, with x 
1 f 1/2 x 3 / 2 x 5 ' 2 hx 7 ' 2 21x 9 / 2 

This series can be plotted against the exact values, and 
in ten thousand even for d as low as five. 



= l/d: 

+ 0{x ll ^)\. (10) 
is accurate to one part 



5 A recursion relation for k(d) 

We can find a simple recursion relation for k{d + 1) in terms of k(d). Observe 
that 

M d M d+1 = f ^ Ty , (11) 
where y = \{d — 1). Using the identity that T(l + y) — yF(y), we have 

M d M d+1 - -Aj-. (12) 

Plugging into the above using equation ([9| to relate k(d) and k(d + 1) to, 
respectively, M(d) and M(d+ 1), we obtain the simple recursion 
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One can easily check that, starting with k(3) — 1/4, this formula allows us 
to reproduce the results of equation ([9| , given in Table 1 . 



6 An explicit formula for k(d) 



Using the recursion (equation ([13|) and anchoring it at fc(2) = 1/tt, we can 
derive an explicit formula for k(d), d > 2. Writing out the values of k(d) for d 
from 3 to 10 allows us to observe that the following formulae hold for d > 2: 

1 (d-2)! 



k(d) = — d odd (14) 



2 (d-1)! 
and 

deven (15) 

7T (a - iy.o 

where subscripts O and E denote, respectively, that the factorial is taken 
only over the odd or even numbers. For instance, \Q\e = 10 ■ 8 • 6 • 4 • 2, and 
9!o = 9 • 7 - 5 • 3 • 1. As this is non-standard notation, I present these formulae 
written in terms of products below as well (for d > 2). 

(d-3)/2 

*w=2 n doM ^ 

n=0 



and 

2n 



n=0 
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